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ABSTRACT 

High-dispersion spectra of 333 stars in the open cluster NGC 6819, obtained using the HYDRA 
spectrograph on the WIYN 3.5m telescope, have been analyzed to determine the abundances of 
iron and other metals from lines in the 400 A region surrounding the Li 6708 A line. Our spectra, 
with signal-to-noise per pixel (SNR) ranging from 60 to 300, span the luminosity range from the 
tip of the red giant branch to a point two magnitudes below the top of the cluster turnoff. We 
derive radial and rotational velocities for all stars, as well as [Fe/H] based on 17 iron lines, [Ca/H], 
[Si/H], and [Ni/H] in the 247 most probable, single members of the cluster. Input T e ff estimates 
for model atmosphere analysis are provided by (B — V) colors merged from several sources, with 
individual reddening corrections applied to each star relative to a cluster mean of E(B — V) = 0.16. 
Extensive use is made of ROBOSPECT, an automatic equivalent width measurement program; its 
effectiveness on large spectroscopic samples is discussed. From the sample of likely single members, 
[Fe/H] = —0.03 ± 0.06, where the error describes the median absolute deviation about the sample 
median value, leading to internal precision for the cluster below 0.01 dex. The final uncertainty 
in the cluster abundance is therefore dominated by external systematics due to the temperature 
scale, surface gravity, and microturbulent velocity, leading to [Fe/H] = —0.02 ± 0.02 for a sub-sample 
restricted to main sequence and turnoff stars. This result is consistent with our recent intermediate- 
band photometric determination of a slightly subsolar abundance for this cluster. [Ca/Fe], [Si/Fe], and 
[Ni/Fe] are determined to be solar within the uncertainties. NGC 6819 has an abundance distribution 
typical of solar metallicity thin disk stars in the solar neighborhood. 

Subject headings: open clusters and associations : individual (NGC 6819), stars : abundances 


1. INTRODUCTION 


NGC 6819 is an old (2.3 Gyr) open cluster whose fun¬ 
damental properties have garnered increasing attention 
in recent years. As detai led in |Anthony-Tw arog, Deliyan-| 
nis, & Twarog (2014|), (hereinafter Paper I), its age 


places it in a sparsely populated range of open cluster 
ages, making it an invaluable testbed for stellar models 
of stars just above the Sun’s mass. It was chosen as a key 
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cluster in our program to map the evolution of Li among 
stars of varying mass as they evolve from the main se¬ 
quence to the tip of the giant branch and beyond. Delin¬ 
eation of such a map requires reliable estimates of stellar 
temperature, luminosity, and metallicity, especially iron, 
and, indirectly, the cluster reddening, distance modulus 
and age. 

Data on individual cluster stars and the surrounding 
field have expanded due to the cluster’s location in the 
Kepler field, bringing the added potential for asteroseis- 
mic insight into the structure of individual stars in a 


variety of evolutionary states (Gilliland et al. 2010). 


Reliable radial velocities are available for a large samp 


of stars e xtending well below the cluster turnoff (Hole et 


al. 2009) (herei nafter H09) and i mprov ed proper-motion 
memberships by Platais et al. (2013) (hereinafter PL) 
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have supplanted the olde r data o f Sanders ( 1972 ). Com¬ 
prehensive broad-band ([Kosvick fe VandcnBerg 1998 


Kalirai et al. 2001; Yang et ai. 20l3f and intermediate- 
band (Paper I) photometric surveys allow precise deter- 


(Paper Ij photometric surveys allow precise 
mination of each star’s position in the color-magnitude 
diagram (CMD), while definitively demonstrating (PL, 
Paper I) that the cluster suffers from variable foreground 
reddening with a range of A E(B — V") = 0.06 mag, 
not surprising given the cluster’s galactic coordinates 
(/, b = 74°, +8°) and a distance of 2.4 kpc from the Sun. 
Until recently, the weakest link in the discussion of the 
properties of the cluster and its stars has been the cluster 
metallicity. As detailed in Paper I, metallicity estimation 
from cluster members has been plagued by small num¬ 
ber statistics often coupled with incorrect assumptions 
of uniform and/or anomalously high reddening, leading 
to a consensus view that NGC 6819 was metal-rich, with 
[Fe/H] ~ +0.1. Precise intermediate-band photometry 
on the extended Stromgren system (Paper I) strongly 
contradicts this claim, producing [Fe/H]= —0.06 ± 0.04 
from 278 single, unevolved F stars. 

Although a fuller discussion of the Lithium abundances 
for our la rge samp le in NGC 6819 is_ in prepa r atio n 
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M5| 7 one preliminary result - the identification of a Li- 

rich giant star - was reported in |Anthony-Twarog et al. 

(2013). We employ a combination of spectral synthesis 
and curve-of-growth analysis to convert measured equiv¬ 
alent widths for the Li line at 6707.8 A to Lithium abun¬ 
dances, A(Li). Our curve-of-growth analysis method de¬ 
pends on subtracting an estimated contribution by an 
iron line at 6707.4 A from the equivalent width, for which 
accurate temperature and iron abundance information is 
required. Given the criticality of a precise estimate of 
the stellar Fe abundance for derivation of Li from spec¬ 
troscopy, especially among the giants, it was decided that 
all measurable lines within the wavelength range of our 
HYDRA high dispersion spectra of the Li 6708 A re¬ 
gion would be used to estimate the cluster metallicity 
for Fe and other elements. In addition to testing the 
photometric abundance, following the pattern laid out in 
our ear lier investi g ations o f the open clusters NGC 3680 
(I Anthony-Twarog et al. 2009) and NGC 6253 (Anthony^ 
Twarog et al. |2010[), such data could reveal the effects of 


stellar evolution on surface abundances through compar¬ 
ison of the red giants to the turnoff stars, while pinpoint¬ 
ing the place of the cluster within the global context of 
galactic chemical evolution. 

The layout of the paper is as follows. Sec. 2 details 
the spectroscopic observations and reduction of images 
to uniformly wavelength-calibrated and continuum-fitted 
spectra. Sec. 3 explains the collation of radial-velocity, 
photometric, proper-motion, and reddening data used to 
define the crucial temperatures, surface gravities, and 
microturbulent velocities for single-star cluster mem¬ 
bers that, together with appropriate model atmospheres, 
translate the spectroscopic measures to abundances. Sec. 
4 supplies a detailed examination of ROBOSPECT, an 
automated line-measurement program used for this study 
which makes internally consistent and rapid abundance 
estimation for hundreds of stars feasible. Sec. 5 presents 
the abundance results for individual stars and the cluster 
as a whole, delineating the impact of the parametric un¬ 


certainties on our final results, and places the data within 
the context of both stellar and galactic evolution. Sec. 6 
summarizes our conclusions. 


2. SPECTROSCOPIC DATA 


Spectroscopic data were obtained for target stars in 
NGC 6819 using the WIYN 3.5-meter telescope [^] and 
HYDRA multi-object spectrograph over 13 nights from 
September and October 2010, June 2011 and February 
2013. The positions of the target stars within the clus¬ 
ter (V,B — V ) CMD are shown in Figure 1. Six con¬ 
figurations were designed to position fibers on a total 
of 333 stars. The brightest configurations include stars 
near V ~ 11, while the fainter configurations reach to 
V ~ 16.5. Individual exposures ranged from 10 to 90 
minutes, with accumulated totals of 2.5 to 4.5 hours for 
stars in the brightest configurations, 10 to 13 hours for 
the configuration with stars of intermediate brightness, 
and over 14 hours total for the faintest configurations. 
Our spectra cover a wavelength range ~ 400 A wide cen¬ 


tered on 6650 A, with dispersion of 0.2 A per pixel. Ex¬ 
amination of Thorium-Argon lamp spectra indicates that 
the line resolution comprises 2.5 pixels, yielding a spec¬ 
tral resolution over 13,000. 

The data were processed using standard reduction rou¬ 
tines in IRAlfj These routines included, in order of ap¬ 
plication, bias subtraction, division by the averaged flat 
field, dispersion correction through interpolation of the 
comparison spectra, throughput correction for individual 
fibers using daytime sky exposures in the same config¬ 
uration, and continuum normalization. After flat field 
division and before the dispersion correction, the long- 
exposure program images were cleaned of cosmic rays 
using “L. A. Cosmic”p (van Dokkum 2001). Real-time 
sky subtraction was accomplished by using the dozens of 
fibers not assigned to stars and exposed to the sky for 
each integration. Composite spectra for each of the six 
configurations were constructed by additive combination. 

The signal-to-noise ratio per pixel (SNR) may be esti¬ 
mated two ways: first, by direct inspection of the spec¬ 
tra within IRAF’s SPLOT utility, using mean values and 
r.m.s. scatter from a relatively line-free region or, second, 
by construction from output files of the ROBOSPECT 
software suite. We found the ROBOSPECT values to 
be entirely consistent with hand-measured SNR values 
and quote the ROBOSPECT-derived values in Table 1, 
computed from the relatively line-free 6680-6694 A re¬ 
gion. SNR estimates reflect the statistics characterizing 
the summed composite spectra. 


3. STELLAR PROPERTIES 

3.1. Radial Velocities and Proper Motions 

An initial spectroscopic sample of probable cluster 
members was taken from the valuable radial-velocity 


3 The WIYN Observatory is a joint facility of the University of 
Wisconsin-Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatory. 

4 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

5 http://www.astro.yale.edu/dokkum/lacosmic/, an IRAF 
script developed by P. van Dokkum (van Dokkum 2001); spec¬ 
troscopic version. 
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survey of NGC 6819 by H09. All stars brighter than 
V ~ 16.75 with radial-velocity membership probabili¬ 
ties greater than 50% were identified as candidates for 
the present study. Stars classed as double-lined spec¬ 
troscopic binaries were eliminated; single-lined systems 
were retained since the existence of the companion would 
have minimal impact on line measurement. Stars were 
not eliminated based upon their position in the CMD 
to avoid biasing the sample against stars undergoing 
anomalous evolution. 

Individual stellar radial velocities were derived from 
each summed composite spectrum utilizing the Fourier- 
transform, cross-correlation facility FXCOR in IRAF. In 
this utility, program stars are compared to stellar tem¬ 
plates of similar effective temperature ( T e ff ) over the 
wavelength range from 6575 A to 6790 A, as well as a 
narrower region in the vicinity of Ha alone. Typical un¬ 
certainties in the individual radial velocities were esti¬ 
mated at 1.15 km/sec. The FXCOR utility also provides 
measures of the line widths within each spectrum, from 
which rotational velocities may be inferred. Observations 
of a pair of radial-velocity standards were obtained dur¬ 
ing each run and processed using the same procedure ap¬ 
plicable to the cluster. Comparison of the velocity zero- 
points using standar d value s from the General Catalog of 
Radial Velocities (Wilson 1953) allowed transformation 
of the cluster data to the standard system, within the 
uncertainties of the measurements. 

As a first check on our spectroscopic data, we can com¬ 
pare our measured radial velocities with those of H09 to 
identify discrepant stars or long-term variables. Elimi¬ 
nating 29 stars classed as spectroscopic binary members 
of NGC 6819, the remaining 304 single stars have a mean 
radial velocity of 2.65 ± 1.36 (s.d.) km/sec. The same 
sample from H09 has a mean velocity of 2.38 ± 0.99 (s.d.) 
km/sec. The dispersion in the residuals between these 
samples is 1.06 km/sec; the predicted dispersion from the 
quoted errors for each star is 1.19 km/sec. The residu¬ 
als for only three stars (WOCS 1007, 2016 and 56018) 
fall more than three sigma from the mean; one of these, 
WOCS 2016, has unusually broad lines (presumably due 
to high rotational velocity) which leads to a larger than 
average uncertainty in the final radial-velocity estimate. 
By contrast, the sample of 29 spectroscopic binaries has 
a mean cluster velocity of 3.45 km/sec and a dispersion of 
9.8 km/sec; the same stars from H09 exhibit a dispersion 
of 9.78 km/sec. We conclude that all single star members 
as classified by H09 are confirmed as such by our data. 
Table 1 contains the derived mean radial and estimated 
rotational velocities for each star in our survey. 

At the start of our program, the only NGC 6819 
prop er-motion study available was that of Sanders”] 
(1972), which proved inadequate for reliable identifica¬ 


tion of probable members, particularly at the fainter 
limit of interest. Fortunately, the comprehensive sur¬ 
vey by PL covers the appropriate range in both area and 
depth, generating probabilities for all stars in our sam¬ 
ple except one. Of 332 stars selected via radial velocity, 
59 have proper-motion membership probabilities below 
50%, 43 of which are in single digits. The PL values for 
individual stars may be found in Table 1. The different 
stellar categories are identified in the CMD, Fig. 1. 


3.2. Effective Temperature, Surface Gravity, and 
Microturbulent Velocity 

In keeping with our approach to spectroscopic abun¬ 
dance determination for previous clusters in this pro¬ 
gram, our default scheme for determining model atmo¬ 
sphere input temperatures is based upon photometric 
color, specifically B — V. Four sources of broad-band 
BV data exist. As a first step, we merged the CCD pho¬ 
tometry presented by H09 as PHOT98 and PHOT03. 
From 923 stars brighter than V = 16.7, excluding 5 stars 
with absolute residuals greater than 0.10 mag, the mean 
offsets, in the sense (PHOT98 - PHOT03), in V and 
B — V are —0.005 ± 0.022 and +0.003 ± 0.022, respec¬ 
tively. While the dispersion in both sets of residuals was 
satisfactory, further analysis revealed that a large frac¬ 
tion of the scatter in V is the result of a nonlinear ra¬ 
dial gradient among the residuals, reaching AV = —0.05 
mag for stars near the core of the cluste r but +0 . 01 for 
st ars in the outer regions of the frame. |Milliman et al.| 
(2014) have independently discovered the same effect 
m a compari son between the VI photometry of |Yang| 


et al. (2013) and PHOT03, implying that the primary 
source of the trend in the current comparison must lie 
with the PHOT03 database. Additional evidence in sup¬ 
port of this conclusion comes fro m a direct comparison 
between the V magnitudes from Yang et al. (2013) 
and the intermediate-band data from Taper 1 where the 
residuals do not exhibit a radial dependence. More im¬ 
portant from our perspective is the absence of a radial 
trend among the residuals in B — V between the two 
samples, PHOT98 and PHOT03. Application of a small 
color term, (B — V ) 0 3 = 0.986(7? — W) 98 + 0.006, reduces 
the scatter in the residuals to ±0.021 mag. The PHOT98 
B — V data were converted to the PHOT03 system and 
averaged for stars common to both samples. 

Next, the composite H09 d ata were tran sformed to 


our adopted standard of Rosvick & VandenBerg (1998). 
From 470 stars with V brighter than 16.5 common to the 
two datasets, removing 10 stars with residuals in B — V 
greater than 0.05 mag, the mean residual, in the sense 
(RV-H09), is +0.006 ± 0.014 mag. A weak color de¬ 
pendence among the residuals was found, (B — V)nv = 
1.007(7? — V)ho 9 , and applied prior to the merger of the 
two databases, eliminating the small zero-point offset. 


Finally, a transformation was derived between the 

B — V indices of 

Kalirai et al. 

(2001 

) and the previ- 

ously merged 7? — 

V data on the J 

kosvic 

k & VandenBerg 


& VandenBerg (fl998j) and |Kalirai et al. 
a ir 


( 200l| ) have 


already been discussed in Paper 1. 

Eliminating 11 stars with absolute residuals greater 
than 0.1 mag and applying a small color term, (7? — 
V)rv = 1.025(7? — V)ka~ 0.018, 846 stars brighter than 
V = 16.5 exhibit a mean residual in B—V of 0.000±0.018 
mag. 

All stars with absolute residuals greater than 0.05 mag 
were individually checked. If 7? — W estimates were avail¬ 
able from more than two sources and one source was 
clearly the origin of the discrepancy, that value was 
dropped. If only two sources of 7? — V existed, an in¬ 
dependent check on the predicted 7? — V was attem pted 
using the published V — I data of Yang et al. (2013). In 


cases where no resolution of the discrepancy was possible 




























































4 


Lee-Brown, Anthony-Twarog, Deliyannis, Rich, Twarog 


the averaged value from all sources was retained. Aver¬ 
aged B — V indices from the four transformed primary 
data sets are given in Table 1 along with the number of 
sources and each star’s estimated individual reddening 
correction. 

As noted earlier, PL identified and mapped variable 
reddening across the field of NGC 6819, a result con¬ 
firmed in Paper I. We used the map of individual red¬ 
dening values derived by PL to estimate by spatial in¬ 
terpolation the degree of reddening affecting each of the 
stars in our spectroscopic sample. Individual reddening 
estimates, with a range of ±0.033 about our adopted 
mean value of E(B — V) = 0.16, were applied to each 
star’s (B — V) color. A temperature for each star based 
on its dereddened (B — V ) color was then derived using 
two primary color-temperature calibrations. 


For dwarfs, we cont inue to us e a calibration (Deliyan- 
n is, Steinhauer fe Jeffries |2002 ) consistent with previous 
spectroscopic studies by this group, namely: 

T e ff = 8575 - 5222.7 (B - V) 0 + 1380.92 (B - V)l + 
701.7 (B - y) 0 [[Fe/H] - 0.15] K. 


For giants, the T e j f -color-[Fe/H] calibration of Ramirez 
& Melendez (2005k was used. We note that for 13 stars 


it was possible to obtain temperatures from both color- 
te mperature c ali brations. Temperatures derived from 


the Deliyannis, Steinhauer & Jeffries (2002) calibration 
are 41 ± 6 K higher tor these mostly subgiant stars. 

Surface gravity estimates (log g ) were obtained by di¬ 
rect comparison of V magnitudes and B — V colors for 
our sa mple of 333 stars to isoch rones from the Y 2 compi¬ 
lation ( |Demarque et al. j|2004 ), constructed for a scaled 
solar composition with [Fe/HJ = -0.06 and an age of 2.3 
Gyr, essentially the same as the comparison presented 
in Paper I. The isochrone’s predicted magnitudes and 
colors were adjusted to match the cluster’s reddening, 
E[B — V) = 0.16 and apparent distance modulus, 12.40. 
For stars that appear to be blue stragglers, surface grav¬ 
ities were estimated by comparing photometric informa¬ 
tion to a grid of younger isochrones of similar composi¬ 
tion. 

We are fortunate to have independent confirmation of 
the log g values for the giant branch stars, thanks to their 
status as Kepler objects of interest. We compared log g 
values from our isochrone match to log g values inferred 
using asteroseismology from Kepler data, as presented 


by Basu et al. (2011) for 21 giants. The seismological 


gravities are m remarkable agreement with the isochrone- 
inferred gravities, differing by an insignificant amount, 
—0.02 ± 0.04 where the listed error is the standard devi¬ 
ation. 

Input estimates for the microturbulent velocity param¬ 
eter were constructed using various prescriptions. For 


formula of Edvardsson et al. 

(1993| was used; a simi- 

far formufation by Ramirez, A 

iende Prieto, & Lambert 


(2013) extends to slightly lower temperatures and grav- 


ities and was used for some stars. For giants, a gravity- 
dependent formula, V t = 2.0 — 0.2 log g, was used. For 
some subgiants and candidate blue stragglers, no suitable 
formula for V t was found other than the purely gravity- 
dependent expression employed by the SDSS collabora¬ 
tion in their DR10 data release and discussion of the 


4. SPECTROSCOPIC PROCESSING WITH 
ROBOSPECT 

With the expansion in spectroscopic sam ples from 
a few dozen stars in previo us cluster work (lAnthony- 


Twarog et al. 2009 2010) to a few hundred in this 


and future analyses, manual measurement of equivalent 
widths (EW) for individual lines is now prohibitive. To 
overcome this obstacle, the automated line-mea suring 
program, ROBOSPECT (|Waters & Hollek ||2013) (here¬ 
inafter WH), has been utilized. For full details, the 
reader is referred to WH; a brief outline of the pro¬ 
gram operation and our procedure is provided here. RO¬ 
BOSPECT measures EW by first determining the con¬ 
tinuum level and constructing a noise profile across the 
wavelength range of interest using an iterative process. 
With the initial continuum and noise levels established, 
spectral lines are tagged at locations specified by a user- 
supplied line list, and other potential lines are automat¬ 
ically identified based on significant deviations of the lo¬ 
cal spectrum from the continuum. For the automatic line 
identification, the significance threshold is user-specified 
and based on the current noise solution. After lines 
are identified though these two processes, they are sub¬ 
tracted from the spectrum and the continuum process is 
repeated and the noise level refined. Through multiple 
iterations, the best-fit continuum and individual line so¬ 
lutions are reached. For this study, each spectrum was in¬ 
dividually corrected in ROBOSPECT for radial velocity 
and run through 25 iterations of continuum fitting and 
line estimation using a gaussian line profile with three- 
sigma automatic line identification and no least-squares 
line deblending. All other parameters for the program 
were set to default values. 

An issue of obvious concern with any automated pro¬ 
cedure is the potential for unreliable EW due to low SNR 
spectra, line blending, or inaccurate radial-velocity cor¬ 
rection. With interactive measurement, such issues can 
be flagged and corrected or eliminated but, with RO¬ 
BOSPECT, flawed EW measures propagate to followup 
steps in the abundance determination process, poten¬ 
tially distorting the final results. It is expected that the 
uncertainty in EW will increase with decreasing SNR, 
potentially reaching a level where the subsequent uncer¬ 
tainty in a star’s calculated abundance renders it useless. 
Similarly, we can expect the EW uncertainty to increase 
as individual lines suffer increased blending from neigh¬ 
boring lines. Finally, a significant error in radial velocity 
will cause ROBOSPECT to misidentify the wavelength 
of line centers. To ensure that our abundances are mini¬ 
mally affected by these issues, a variety of tests were con¬ 
ducted to determine if our line list and ROBOSPECT 
input parameters lead to robust results over the range 
in SNR and radial-velocity uncertainties spanned by our 
spectra. Since ROBOSPECT will be used in future anal¬ 
yses with similar input parameters, the level of detail 
supplied is somewhat greater than usual. 

The line list used during this study was constructed by 
visually identifying relatively isolated, medium-strength 
(EW ~ 10 to 200 mA) li nes in the so lar spect rum, using 
the solar atlas of Wallace, Hinkle, & Livingston (1998) 
as a guide. Atomic information for each fine was re¬ 


abundance analysis pipeline for APOGEE spectra (Ahn 
et al. |[2014| . 


trieved from the VALD database (Kupka et al. 1999). 


Log gf values given by VALD were then modified to 
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force-fit abundances from ROBOSPECT analysis of so¬ 
lar spectrum EW to abundances in the 201 0 version of 
the abundance analysis software, MOOG (Sneden 1973). 
We used day-time sky spectra, taken as part of the caii- 
bration data sets, as solar spectra for this normalization. 
MOOG analysis w as cond ucted using the abfind driver 
and a solar model (Kurucz 1995) with the following at¬ 
mospheric parameters: ('1 ' e ff , log g,v t , [Fe/H ]) = (5770 
K, 4.40, 1.14 km/s, 0.00). 

Our final line list contains 22 lines of interest (17 Fe, 3 
Ni, 1 Ca, 1 Si), presented along with the relevant atomic 
parameters in Table 2. To minimize the impact of blend¬ 


ing, all selected lines are at least 0.5 A from any line with 
EW greater than 5 mA in the solar spectrum. To further 
minimize measurement distortions in spectra with lower 
SNR, the line list input for ROBOSPECT included any 
significant feature within 2 A of a line of interest. EW 
for these features were disregarded during the abundance 
analysis. The extraneous entries were added after observ¬ 
ing unusual amounts of scatter among the EW for lines 
with close neighbors, particularly at lower SNR. 

To evaluate the performance of our line list and to 
quantify the effects of a spectrum’s SNR on EW values, 
we tested ROBOSPECT on 100 solar spectra, 25 each 
with mean SNR of 160, 130, 95 and 70. These represent 
typical SNR for our program stars (see Fig. 2). We note 
that WH include an evaluation of ROBOSPECT’s per¬ 
formance over a wide range of line strengths and SNR; 
our test results are consistent with theirs. The mean 
standard deviations in EW for our four test samples were 


4.2, 4.7, 7.5, and 11.1 mA, in order of decreasing mean 
SNR. The uncertainty this scatter introduces into our 
abundance results is small relative to the uncertainty due 
to external systematics. At all SNRs tested, we observed 
no significant correlation between the mean EW of a line 
and its associated uncertainty. Taken together, these re¬ 
sults indicate that our line list is robust over the range 
in mean EW and SNR of interest. 

The effect of potential errors in the determined radial 
velocity on EW measurement was evaluated by artifi¬ 
cially increasing and decreasing our best estimate of ra¬ 
dial velocity by a flat amount in the set of 100 solar spec¬ 
tra. For a shift of ±2.5 km/s, corresponding to a radial- 
velocity uncertainty greater than that of 99% of our pro¬ 
gram stars, the mean deviation in EW between the true 
and adjusted spectra was —0.7 mA. When the shift was 
increased to 5.0 km/s, corresponding to a radial-velocity 
uncertainty greater than that of all program stars (in¬ 
cluding the rapid rotator WOCS 2016), the mean devi¬ 
ation in EW increased to —3.5 mA. In both cases, RO¬ 
BOSPECT correctly identified essentially the same num¬ 
ber of entries from the line list as it had with the correct 
radial velocity. We conclude that the effect of radial- 
velocity error on our abundance results is negligible. 

Finally, we compared ROBOSPECT results with man¬ 
ual EW measurements for 18 random program stars us¬ 
ing a 10-line subset of our line list. Manually obtained 
EWs were measured using IRAF’s SPLOT tool. RO¬ 
BOSPECT’s results largely agree with the manual val¬ 
ues, with ROBOSPECT’s EW on average being 7.7 ±4.2 
mA (s.d.) lower than the manual results. We note that 
this consistent underestimate of EW by ROBOSPECT 
(or, over-estimate using SPLOT and interactive measur¬ 


ing) is similar in solar spectra and in the spectra of pro¬ 
gram stars, and is not a function of SNR to any de¬ 
tectable extent, so we anticipate no effect on our abun¬ 
dance estimates based on differential analysis with re¬ 
spect to the sun. 

Using our measured radial velocities and constructed 
line list, EWs were measured for all 333 stars in our 
observing program. We discarded negative EW (RO¬ 
BOSPECT’s designation of emission lines) due to non- 
convergent fitting solutions, artifacts of cosmic ray re¬ 
moval, large noise spikes, or nonexistent lines in the mea¬ 
sured spectrum. Approximately 4% of the measured lines 
of interest returned negative values. 


5. ABUNDANCE DETERMINATIONS 

Model atmosphere s were con structed for each program 
star using the grid of Kurucz (1995) and the input T e //, 
log g and microturbulent velocity values listed in Table 
1. As a reminder, each star’s T e ff estimate is based 
on its B — V color, individually dereddened based on 
its spatial position within the field of the cluster. For 
any star that is not actually a cluster member, this spa¬ 
tially interpolated reddening correction may be spurious; 
the gravity estimate based on an isochrone matched to 
the cluster CMD almost certainly would be as well. For 
suspected binaries, the color might reflect the combined 
light of the two stars and affect the assigned temperature. 
While this might not be a large issue for SB1, a few can¬ 
didate SB2 were highlighted by examination of spectra: 
WOCS 7009, 26007, 35025, 9004 and 60021. As a conser¬ 
vative approach, we restricted our analysis to stars with 
membership probablity > 50% (PL) and no evidence of 
binarity. We also excluded from our final analysis several 
stars with (B — U)o greater than 1.35, for which severe 
line-blending makes reliable measurement of equivalent 
widths problematic. This left us with a sample of 247 
single, member stars. 

Each star’s measured equivalent widths and model 
serve as input to the abfind routine of MOOG to pro¬ 
duce individual [A/H] estimates for each measured line 
for each star. A large volume of spectroscopic informa¬ 
tion is available for interpretation from over 7000 mea¬ 
sured lines, with up to 17 [Fe/H] values for each star. 
Before constructing median [Fe/H] values for each star 
and for the sample as a whole, the following filters were 
applied to the individual abundance estimates: abun¬ 
dances based on equivalent width measurements smaller 
than three times the expected error in equivalent width 
for each star based on the star’s SNR, were not con¬ 
sidered, nor were equivalent widths for very strong lines 
(EW > 200 mA). We note that the consequences of using 
200 or 250 mA as the upper limit for EW are negligible. 
Estimated abundances from individual lines deviant by 
more than 1 dex from solar were not included in the de¬ 
termination of [Fe/H] for each star. 

Figures 2 and 3 show the spread of stellar [Fe/H] esti¬ 
mates for each star as a function of the individual spec¬ 
trum SNR and unreddened B — V color for each star. 
A few outliers with SNR below 50 are obvious, as are a 
few of the brightest (and coolest) stars with higher than 
typical abundance. Figure 3 suggests that abundances 
for the cooler stars appear to be slightly below the range 
for the stars nearer the turnoff. We further divided our 
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sample of 247 likely members into stars bluer than and 
redder than (B — V)o = 0.60 to separate out the cooler 
and more evolved stars and re-examined the effect of fil¬ 
ters that precede the construction of each star’s median 
[Fe/H] abundance. Do some of our 17 iron lines produce 
abundances in dwarfs that are unacceptably noisy? 

By examining the effect of imposing different criteria 
for excluding small EW measures, we were able to high¬ 
light several iron lines, which for dwarfs, lead to biased 
[Fe/H] values in the sense that the abundance values for 
these lines in dwarf stars are higher than the ensemble 
median value by several tenths. The reason for this is 
clear; if the median EW for a particular absorption line 
is small, imposing a 3-sigma cut will preferentially pre¬ 
serve stars in the sample with larger EW, and conse¬ 
quently larger [Fe/H] values. We excluded lines in dwarfs 
for which imposing 3-sigma cuts changed the median EW 
for the line by more than 25% in comparison with a more 
generous 1-sigma cut. As a reminder, these lines are mea¬ 
sured in the solar spectrum at considerably higher SNR 
so we may still place considerable confidence in the log 
gf values determined from those lines. 

As this still leaves information from several thousand 
separate iron line measurements, we were able to further 
explore the consequences of these filters. A diagnostic 
diagram is presented in Figure 4, demonstrating the lack 
of a trend in the median [Fe/H] for each of our 17 iron 
lines with wavelength, where the median value is con¬ 
structed from the full sample of 247 probable-member, 
single stars for lines denoted by black symbols; the six 
lines denoted by red symbols represent measures for gi¬ 
ant stars only. The upper panel of this figure shows the 
MAD (median absolute deviation) statistic for each line, 
a robust estimator of variance, constructed by compar¬ 
ing the [Fe/H] value for each star to the median value for 
that particular absorption line. 

For the subset of 247 single member stars bluer than 
(B — V )o < 1.35, the median [Fe/H] value is — 0.03±0.06, 
with abundances for the other elements that are entirely 
consistent with solar: —0.01 ± 0.10,-0.01 ± 0.06 and 
0.00±0.06 for [Ca/H], [Si/H] and [Ni/H], respectively. Of 
these 247 stars, 200 are main sequence and turnoff stars, 
with a median iron abundance of —0.02 ± 0.05, while the 
remaining 47 cool and evolved stars show a lower [Fe/H], 
—0.09 ± 0.05. We note that this apparent discrepancy 
between giant and dwarf [Fe/H] val ues is nearly identi- 
cal to that found in NGC 3680 (Anthony-Twarog et al. 
2009). 

Fot all of these quantities, the indicated error is the 
MAD statistic; for a large sample size drawn from a 
normally distributed population, the normalized MAD 
statistic, MADN = 1.48 MAD, approximates well the 
sample standard deviation, permitting an estimate of the 
traditional standard error of the mean through division 
by (N — l) 1 / 2 . For the sample of 247 single, probable- 
member stars, the MAD statistic implies a value for 
MADN = 0.09; as this statistic is a good proxy for a 
standard deviation, a standard error of the mean of 0.01 
is implied. These values can be confirmed by computing 
averages, standard deviations and s.e.m. values in linear 
space; for these stars, the average [Fe/H], computed in 
linear space and then converted to logarithmic values, is 
— 0 . 02 ^ 0)53 (standard deviations). 


Table 3 includes the median [Fe/H], along with the 
number of iron lines included in the estimate, as well 
as estimates of [Ca/H], [Si/H] and [Ni/H] for 333 stars. 
Stars not part of the set of 247 single, probable mem¬ 
bers stars are flagged with a note indicating non-member 
of binary status. The statistic that accompanies each 
[Fe/H] estimate is the MAD. We note that the 29 stars 
not included in these estimates of cluster [Fe/H] due 
to evidence of binarity, yield a similar median [Fe/H] 
value if the suspected double-lined spectroscopic bina¬ 
ries are excluded. For member stars which are mostly 
SB1, [Fe/H]= —0.05 ± 0.06, consistent with the more 
conservative sample that excludes binaries. 

Conversions of photometric parameters to T e // and 
V t are fundamentally discontinuous for dwarf and giant 
stars; to assess the effects of parameter choices on the 
derived abundances, each star was reanalyzed with at¬ 
mospheric parameters incremented and decremented by 
the following amounts: 100 K for T e ff, 0.25 for log g and 
0.25 km/sec for V). A summary of the extensive results 
is shown in Figure 5 where once again median values 
and variances indicated by median absolute deviations 
are compiled and shown for stars in different color bins. 
This figure illustrates some well-known dependences of 
abundance determinations on stellar parameters, namely 
that main sequence star abundances are very sensitive to 
T e ff estimates, while evolved stars demonstrate a higher 
sensitivity to surface gravity and consequently to micro- 
turbulent velocity estimation schemes. 

Although the discrepancy between the giant and near- 
turnoff samples is small, it is of value to explore modest 
parameter changes that could eliminate it. Reference 
to Figure 5 suggests that errors in log g have relatively 
little impact; we remind the reader that the log g val¬ 
ues employed are consistent with values suggested by 
asteroseismic analyses of NGC 6819 giants. The rela¬ 
tion used for microturbulent velocities for the evolved 
stars is more tenuous; a fairly modest decrement of 0.2 
km/sec would raise the giant abundances to match the 
[Fe/H] representative of the near-turnoff stars. Finally, 
increased abundances for both regions of the CMD would 
accompany higher temperatures, though a simple incre¬ 
ment to the mean cluster reddening will have differing 
effects for the two classes of stars, both because of the 
slight color-dependence of the effect on the abundance as 
shown in Figure 5 and the larger increment implied for 
warmer stars resulting from the same higher reddening. 

To define the final estimate of [Fe/H] for the cluster, 
we default to the dwarfs, where the sensitivity to the 
changes in the three key input parameters (T e ff, Vt, log 
g) is collectively minimized, the color range is small, and 
the parametric differences between the program stars and 
the solar reference spectra are significantly smaller than 
for the evolved stars. For log g and V t we adopt 0.1 dex 
and 0.2 km/sec as generous estimates for the uncertain¬ 
ties. From Fig. 5, these translate into 0.003 and 0.008 
dex, respectively. As usual, the more sensitive parameter 
is the temperature, tied to the B — V color transforma¬ 
tion, which has two primary components, the uncertainty 
in the mean cluster reddening, E(B — V), and the uncer¬ 
tainty in the photometric zero point. From Paper I, the 
photometric internal and external errors combined lead 
to an uncertainty of 0.007 in E[B — V), while the pho¬ 
tometric merger of the four broad-band surveys above 
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implies an uncertainty in the B — V zero-point of 0.004 
mag. Combined, the (B — V) 0 scale is reliable to 0.008 
mag. On our temperature scale, this translates to ap¬ 
proximately 32 K or, from Fig. 5, 0.012 dex in [Fe/H], 
Combining all three parameters, one arrives at a final 
external uncertainty of 0.015 dex, much larger than the 
internal precision of the spectroscopic average. From 200 
turnoff stars, the mean [Fe/H] for NGC 6819 is deter¬ 
mined to be —0.02 ± 0.02. 

Among the prior determinations of NGC 6819’s heavy 
element ab undance, the result published by Bragaglia ct 


al. (2001) has anchored a string of moderately super- 
solar values. Their study of three clump giants uti¬ 
lized the high dispersion (R ~ 40, 000) spectrograph on 
the Galileo Italian National Telescope (TNG). Though 
the sample size was small, their analysis incorporated 
hundreds of lines, including nearly 100 iron lines, six 
of which were Fe II lines to permit an accurate ver¬ 
ification of surface gravity. The reddening value of 
E (B — V) = 0.14 they adopted, however, was some¬ 
what smaller than the mean value derived in Paper I 
and utilized here, E (B — V) = 0.16. We repeated pa ram¬ 
eter determinations for the three giants studied by |Bra-| 


gaglia et al. (2001), including independent determina¬ 
tions of spatially-dependent reddening and T e // values, 
log g values from isochrone comparisons and microturbu- 
lent velocity values based on the surface gravities. Using 
the estimates of A [Fe/H] for each of the principal at - 
mospheric parameters cited by Bragaglia et al. (2001), 
their equivalent width analysis would yield an average 
abundance 0.07 dex lower with our atmospheric param¬ 
eters, lowering their [Fe/H] for the cluster from +0.09 to 
+0.02 ± 0.03, essentially solar. 


6. DISCUSSION AND CONCLUSIONS 

High dispersion spectra of 333 potential members of 
the old open cluster, NGC 6819, have been processed 
and analyzed. From the wavelength region near Li 6708 
A, abundances have been obtained for Fe, Ni, Ca, and 
Si for 247 highly probable, single-star members ranging 
from the tip of the giant branch to well below the top 
of the CMD turnoff. Using up to 17 Fe lines per star, 
the cluster exhibits a slightly subsolar [Fe/H], a conclu¬ 
sion which remains unchanged if the data are sorted into 
evolved stars (47) or turnoff/unevolved stars (200). Due 
to the cumulative impact of a few thousand Fe lines, the 
standard error of the mean for [Fe/H] tied to internal 
errors is below 0.01 dex, leaving the systematics of the 
temperature, log g , and microturbulent velocity scales 
as the overwhelming source of uncertainty in the final 
cluster metallicity. Once again, thanks to the size of the 


sample, we can minimize the impact of these parame¬ 
ters by looking only at the 200 unevolved stars since the 
red giant abundances have a higher degree of sensitiv¬ 
ity to microturbulent velocity. From these stars alone, 
[Fe/H] = —0.02 ± 0.02, where the uncertainty includes 
both internal and external errors. The lower metallicity 
compa red to the commonly a dopted value of [Fe/H] = 
+0.09 (Braga glia et al. (2001 ) continues a trend first in¬ 
dicatedDythelntermematfrUand photometry of Paper I, 
where precision photometry was used to derive a reliable 
cluster reddening and confirm the existence of variable 
reddening across the face of the cluster (PL), leading to 
subsolar [Fe/H] from either m\ or hk indices. 

Within the context of Galactic evolution, NGC 6819 
exhibits no features which distinguish it from a typical 
thin-disk population formed within the last 5 Gyrs within 
1 kpc of the Sun’s galactocentric distance. Its [Fe/H] 
places it at the mean for open clusters at the solar galac¬ 
tocentric radius, a distribution which exhibits a scatter 


at on ly the ±0.1 dex level (Twarog, Ashman, & Anthony- 


Twarog 1997). Its abundance ratios relative to Fe are all 


consistent with solar, a pattern which holds true within 
the scatter for Ni, Ca, and Si for stars with [Fe/H] ^0.0, 
whether they are cla ssed as members of t he thin or thick 
disk population (see Hinkel et al. (2014) and references 
therein). NGC 6819 forms a rich ancTideal link for tests 
of ste llar evolution theory midway in age between NG C 
5822 (Carraro et al. |20lip and M67 ( |Sandquist |2004[ , a 
characteristic which will prove valuable in the analysis of 
the evolution of Li for intermediate mass stars evolving 
off the main sequence and up the giant branch. 
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Fig. 1. — Color-magnitude diagram of spectroscopic sample. Blue open circles denote single, probable member stars. The few large 
magenta asterisks mark stars with very cool colors, while black crosses mark probable nonmembers and red filled squares denote probable 
member spectroscopic binaries. 
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Fig. 3.— For 247 single, probable member stars in NGC 6819, median iron abundances relative to solar is illustrated as a function of 
reddening-corrected color, (B — V)o- 
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Fig. 4. — For single, probable member stars, the median abundance is shown for each of the 17 iron lines used, plotted as a function 
of wavelength expressed in A. The upper panel shows the median absolute deviation of values about each line’s median value. Values 
denoted by red symbols indicate lines for which only giant stars contribute to the median abundance value. In the bottom panel, the 
dashed horizontal line denotes the cluster [Fe/H] value of -0.03. 
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Fig. 5.— Median values of A[Fe/H] for stars in unreddened color bins 0.2 mag wide illustrate the effects of increasing the model atmosphere 
temperature by 100 K (blue filled squares), of increasing the microturbulent velocity parameter by 0.25 km/sec (red cross symbols) and 
lowering the surface gravity parameter, log g by 0.25 (green open triangles). 
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Table 1. Stellar Parameters for NGC 6819 Sample 
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Table 2. Lines used for Abundance Analysis 


Wavelength 

Element 

EP 

log gf 

6597.56 

Fe 

4.80 

-1.04 

6609.11 

Fe 

2.56 

-2.71 

6609.68 

Fe 

0.99 

-5.24 

6627.54 

Fe 

4.55 

-1.57 

6646.93 

Fe 

2.61 

-3.87 

6653.91 

Fe 

4.15 

-2.44 

6677.99 

Fe 

2.56 

-1.50 

6703.57 

Fe 

2.76 

-3.02 

6710.32 

Fe 

1.49 

-4.77 

6725.36 

Fe 

4.10 

-2.35 

6726.67 

Fe 

4.61 

-1.16 

6733.15 

Fe 

4.64 

-1.62 

6750.15 

Fe 

2.42 

-2.80 

6806.86 

Fe 

2.73 

-3.14 

6810.27 

Fe 

4.61 

-1.08 

6820.37 

Fe 

4.64 

-1.20 

6837.01 

Fe 

4.59 

-1.82 

6717.68 

Ca 

2.71 

-0.22 

6721.85 

Si 

5.86 

-1.09 

6643.63 

Ni 

1.68 

-1.83 

6767.77 

Ni 

1.83 

-2.23 

6772.31 

Ni 

3.66 

-1.00 


Note. — Wavelengths are expressed in 
A units, excitation potentials in electron 
volts. All abundances refer to neutral 
atoms. 
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Table 3. Abundances for Individual Stars in NGC 6819 Sample 


WOCS ID 

[Fe/H] 

MAD 

N(Fe) 

[Ca/H] 

[Si/H] 

[Ni/H] 

Note 

1002 

-0.31 

0.22 

15 


0.17 

-0.34 

NM/BC 

1004 

-0.20 

0.09 

16 


0.05 

-0.06 


1007 

-0.21 

0.13 

16 


0.29 

-0.03 

NM/BC 

1014 

-0.28 

0.20 

16 


0.10 

-0.10 

NM/BC 

1016 

-0.75 

0.61 

16 

-0.67 


-1.47 

NM/BC 

2001 

-0.07 

0.12 

16 


0.12 

-0.02 

NM/BC 

2003 

-0.13 

0.12 

16 

0.53 

-0.05 

-0.14 


2004 

-0.04 

0.10 

16 


0.23 

0.02 


2006 

-0.11 

0.05 

16 

0.20 

0.05 

0.04 


2007 

-0.07 

0.11 

15 

0.15 

0.07 

0.11 


2012 

-0.13 

0.07 

16 

0.19 

0.06 

-0.04 

NM/BC 

2016 

-0.30 

1.04 

6 



-1.58 

NM/BC 

3001 

-0.16 

0.09 

16 

0.18 

-0.04 

-0.19 


3003 

-0.37 

0.10 

15 

-1.99 

-0.21 

-0.34 


3004 

-0.07 

0.09 

16 


0.12 

0.01 


3005 

-0.16 

0.07 

16 

0.14 

0.03 

-0.07 


3007 

-0.13 

0.08 

15 

0.22 

-0.01 

0.02 


3009 

-0.13 

0.12 

16 

0.26 

0.09 

0.02 


3021 

0.00 

0.11 

14 

0.33 

0.03 

0.02 

NM/BC 

4001 

-0.08 

0.07 

15 

0.15 

0.13 

0.05 




